Phase Space Dynamics of Non-Gravitational Interactions between Dark Matter and 

Dark Energy: The Case of Ghost Dark Energy 
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We study the phase space asymptotics of the so called Veneziano ghost dark energy models. Mod- 
els where the ghost field's energy density: i) pghost oc H, and ii) p g host cc H + H 2 , axe investigated. 
Both, cases with and without additional non-gravitational interaction between cold dark matter and 
ghost dark energy, are subject to scrutiny. We pay special attention to the choice of phase space 
variables leading to bounded and compact phase space so that no critical point of physical interest 
is missing. A rich asymptotic structure is revealed: depending on the kind of non-minimal coupling 
critical points associated with radiation dominance, matter dominance, cold dark matter/ghost dark 
energy scaling, and ghost dark energy dominance, are found. Past and future attractors, as well as 
saddle equilibrium points, are identified in the corresponding phase spaces. 
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I. INTRODUCTION 

The origin and nature of the present stage of accel- 
erated expansion of the universe persists as one of the 
unexplained mysteries in theoretical physics. Assuming 
Einstein's general relativity (GR) as the classical theory 
of the gravitational interactions leads to a large gallery of 
models which explain the current inflationary stage of the 
cosmic evolution within appropriate data accuracy (see 
[l[ for well-known reviews on the sub ject ). The so called 
'ghost' dark energy (GDE) models [3-[l3j have joined the 
last this vast gallery. 1 According to this proposal the 
cosmological constant arises from the contribution of the 
so called Veneziano ghost fields, which are required for 
the resolution of the U(l) problem in the low-energy ef- 
fective theory of quantum chromo-dynamics (QCD) [HI]. 
Although the Veneziano ghosts are unphysical and make 
no contribution in flat Minkowski spacetime, if curved or 
time-dependent spacetime backgrounds are considered, 
the cancellation of their contribution to the vacuum en- 
ergy leave a small energy density p ~ Aq CD H, where 
H is the Hubble parameter and A qcp ~ lOOMeV is 
the QCD mass scale @-IH (see also \m flTjj ) . If assume 
H ~ 10 _33 eV the vacuum energy contribution from the 
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1 From the cosmological point of view ghost dark energy is a par- 
ticular case of dark energy models with generalized equations of 
state [3. 



ghost fields gives the right magnitude ~ (3 x 10~ 3 ey) 4 
[J|. This lucky circumstance is what have stimulated 
the interest of cosmologists in these models of dark en- 
ergy. Their most attractive feature resides in the fact 
that nothing behind the standard model (SM) of parti- 
cles and general relativity is required to explain the origin 
of the dark energy. 

The basic ideas of the model were put forth in @, H[ 
by means of a simplified model. In reference [|| a more 
realistic model was considered and it was shown that 
indeed QCD ghosts produce dark energy in the right 
amount. The cosmological dynamics of a simple GDE 
model where the energy density of the DE is propor- 
tional to the Hubble parameter p g d e — olH (a ~ ^qcd) 
has been studied in Ref.Q. In this model the universe 
has a de Sitter phase at late time and begins to accelerate 
at redshift around z„ ~ 0.6. The authors fit the model 
and give the constraints on the model's parameters with 
current observational data including Snla, BAO, CMB, 
BBN and Hubble parameter data. The impact of ad- 
ditional non-gravitational interaction between GDE and 
cold dark matter (CDM) on the cosmic dynamics was 
also investigated in [5|. It was argued therein that, since 
the squared sound speed of GDE c 2 is negative, this may 
cause an instability. The issue was investigated in detail 
in Ref.|8]. It was found that, due to non-positivity of the 
squared sound speed, both non-interacting and interact- 
ing GDE models are classically unstable against pertur- 
bations in flat and non-flat Friedmann-Robertson- Walker 
(FRW) backgrounds. The seeming instability problem of 
GDE models raised by negativity of c 2 is overcame by 
means of a simple yet solid argument @: unlike other 
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ghost dark energy models, where the ghost field becomes 
a real propagating physical degree of freedom subject to 
stringent constraints (see Ref.|l8[), the Veneziano ghost 
is not a new physical propagating degree of freedom and 
the corresponding GDE model does not violate unitarity, 
causality, or gauge invariance as clearly stated in Q . 

Models of GDE where p gde = aH + /3H 2 (the term 
oc H 2 is inspired by a decaying cosmological constant 
[l9l 1201 ]) were explored in [f| (see also @j). Other as- 
pects of GDE models such as: i) equivalence with kinetic 
k-essence [9j, ii) connection with Brans-Dicke, and f(R) 
theories, and tachyons [T3| , and iii) thermodynamical is- 
sues have been also investigated. In the latter case 
the GDE energy density p gde is shown to be related with 
the radius of the trapping horizon Ft fllj : 

Pgde = —~ = a 1 - e)\ H 2 + -j, 

where e = rx/ZHrT- If neglect spatial curvature, as we 
do in this paper, the trapping horizon is coincident with 
the Hubble horizon Tt = and p gde = oe(l — e)H. 

Here we also neglect the contribution oc e as it is done 
in most works on GDE 

@-d0-[ll- This Particular case, 
strictly speaking, corresponds to dark energy dominance 
so that, in general, the contribution of e to the dynam- 
ics has to be taken under consideration, as it is done 
in In that reference the global dynamical behav- 

ior of the univers e accelerate d by the QCD GDE, where 
Pgde = a(l — e)y ' H 2 + k/a 2 , was investigated by using 
the dynamical systems tools. Depending on the values of 
free parameters only two critical points of physical rel- 
evance were found. 2 Due to the choice of phase space 
variables in [l3j]: p — Q dm/^gde = Pi/3H 2 ) and e, 
both points are correlated with CDM/GDE-scaling be- 
havior since in both cases p = p c 7^ 0. 

The dynamical systems tools supply powerful means to 
extract very useful information out of cosmological mod- 
els without worrying about exact solutions. There is an 
isomorphism between particular solutions of the cosmo- 
logical field equations and points of the equivalent phase 
space. Critical (also equilibrium or fixed) points in the 
phase space can be correlated with particular (classes 
of) solutions of the cosmological equations. If the criti- 
cal point were either a future or a past attractor, then, 
the corresponding cosmological solution were a generic 
or preferred solution of the cosmological equations since, 
quite independent on the initial conditions, orbits in the 
phase space start at a past attractor and end up at a 
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signals bifurcation which is a very interesting finding in Ref . fl3[ . 



future attractor whenever these attractors exist. How- 
ever, no matter how useful these tools are, one have to 
take with care the results of the linear dynamical systems 
study and never should overestimate their reach. It is 
central for the linear study, in particular, that the equiv- 
alent phase space be bounded and compact whenever it is 
possible. Only in this case all of the relevant equilibrium 
points are 'visible'. Otherwise several critical points may 
be missing and one will need of additional phase space 
variables to cover the region not covered by the original 
variables (see section III below). This is precisely what 
we think happens in the study of reference |13j . Actu- 
ally, due to the definition of the variable p (see above), 
the phase space in the mentioned reference is the un- 
bounded plane region {(p, e) 1 < p < 00, e < 1}. Hence, 
several critical points of physical relevance, in particular 
those associated with cosmological solutions which do not 
represent scaling behavior, may be missing in (l3| . 

The aim of the present paper is the detailed investi- 
gation of the phase space dynamics of GDE cosmologi- 
cal models with/without additional non-gravitational in- 
teraction with CDM, by means of the dynamical sys- 
tems tools. We shall consider GDE models where: i) 
Pgde = aH 0, i i£Dl, and ii) p gde = aH + pH 2 $, Q . 
Due to the amount of results to discuss, the most realis- 
tic case when p gde = cx(l — e)H will be left for a sequel 
to this work. Our goal here will be to uncover the main 
features of the asymptotics of these cosmological models 
within the equivalent phase space so that we can corre- 
late equilibrium points (past/future attractors and saddle 
points) with their generic cosmological solutions. We pay 
special attention to the choice of phase space variables 
which lead to bounded and compact physically meaning- 
ful phase space, so that no critical point of cosmological 
interest is missing. The arising asymptotic structure of 
the GDE models is richer than the one found in Ref. [l3[ . 
In the interacting case, depending on the kind of assumed 
interaction, besides the already mentioned CDM/GDE 
scaling equilibrium points, also critical points associated 
with radiation dominance, matter dominance, and GDE 
dominance are found. 

We start by illustrating the method through investi- 
gating the dynamics of a well-known physical system in 
section II: a quintessence model with an exponential self- 
interaction potential. In this context we shall look at 
the role non-gravitational interactions play in the cosmic 
dynamics by comparing with well-known results for the 
non-interacting case [2l| . Non-gravitational interactions 
have proven useful to overcome the coincidence problem 
in the context of GR models of dark energy [22h24| . How- 
ever, the success of the interaction in solving the prob- 
lem depends on the assumed kind of interaction. Here 
we shall consider specific kinds of interaction commonly 
studied in the bibliography (2514271 ]. In section III we 
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give brief tips on (linear) phase space analysis which will 
be useful latter on when we explore the phase space dy- 
namics of GDE models. Then we apply the methodology 
used in section II to explore the ghost dark energy mod- 
els pgde (x H in section IV, and p g d e <x H + H 2 in section 
V. The results of our study and their implications will 
be discussed in section VI, while brief conclusions will be 
supplied in the final section VII. Through the paper we 
use natural units (8nG = Sir/rn 2 PL = c = 1). 



II. INTERACTING QUINTESSENCE 

In this section we shall discuss in a more or less ped- 
agogical way, how the quintessential cosmological dy- 
namics driven by an exponential self-interaction poten- 
tial looks like in the equivalent phase space. Although 
the phase space dynamics of non-interacting quintessence 
fields has been investigated in detail [21| , the phase space 
description of the interacting case is less known (see, how- 
ever [22 - l24| ). In the last instance, the discussion in the 
present section will serve as an illustration of the method 
and terminology we shall apply to the investigation of the 
ghost dark energy case in the next sections. Throughout 
the paper we shall consider Friedmann- Robertson- Walker 
(FRW) spacetimes with flat spatial sections, 

d,s 2 = -dt 2 + a 2 (t)S ik dx l dx k , i,k = l, 2, 3. 

For simplicity we shall assume that the cosmic dynam- 
ics is fueled by a mixture of: i) pressureless CDM with en- 
ergy density p m , and, ii) dark energy. In this section the 
dark energy component is in the form of a quintessence 
field tp, with energy density p v , and parametric pressure 
p v given by: 

Pv = Cp 2 l<l + V{<p), Pv = <j> 2 /2 - V{tp), 

respectively. In these equations V — V(ip) is the self- 
interacting potential for quintessence. We also assume 
that the CDM and the quintessence field interact through 
an additional force of non-gravitational origin. 

The quintessence field equations are the following: 



3H =p m + P v , 
p v + 3H(u 9 + l)p v 

p m + 3Hp rn = Q, 



-Q, 



(1) 



where Q stands for the non-gravitational interaction 
term, and ui v - the quintessence equation of state (EOS) 
parameter: 



_ <p 2 -2V 
p v ~ <p 2 + 2V 



(2) 



Our goal will be to study the dynamics dictated by 
equations (TTJ) in the equivalent phase space. In order 
to obtain an autonomous system of ordinary differential 
equations (ODE) out of ([lj we introduce the following 
dimensionless phase space variables (2lj : 



P 



(3) 



which amount to the (square of the) dimensionless kinetic 
and potential energy density of the quintessence field re- 
spectively. In terms of these variables the Friedmann 
constraint (first equation in JTJ) can be written as 



3H 2 



(4) 



while 



IT 



if 



where the tilde accounts for derivative with respect to 
the parameter, r = J da/ a => dr = din a = Hdt. 

Since < J7 m < 1, then, it follows that, < x 2 + 
y 2 < 1. We shall be concerned here with cosmological 
expansion (H > 0), so that, only non- negative values of 
the variable y (y > 0) will be considered. 

Other physical parameters of interest are: 



Pip 

3H 2 



!f 



Me-**)- 



(5) 



The following autonomous system of ODE is obtained 
out of the original set of cosmological equations: 



x(l + y 2 - x 2 ) 



3 V J± 2 _ Q 
2 V V 6H 3 x' 



-y[l + x 2 -y J 



2V^ 
3 V 



(6) 



Depending on the assumed kind of interaction Q, and 
of the self-interacting potential V, one might need to 
introduce one or two additional phase space variables. 
For convenience the interaction term can be written as 
Q = HS, where S is an energy density ascribed to the 
interaction. If S > 0, the transfer of energy is from the 
DE to the CDM component and vice versa. 

Here we shall concentrate on the exponential potential, 
»-w (y anc | ^ are constants), which is one of 
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the most common potentials found in the applications to 
scalar field models of dark energy. In this case V v /V = 
— /x. We shall assume different forms of the interaction 
which are frequently found in the bibliography on the 
subject [II EE III- 



A. S = So = const. [2 

Here we shall assume positive Sq > 0, so that the trans- 
fer of energy is from the quintessence field into the CDM. 
In this case it is convenient to introduce an additional 
variable, 



3H 2 



3ff 2 + «V 



< z < 1, 



to the variables in ([3]). Notice that the value z = 1 
corresponds to either, i) no interaction (6q = 0), or ii) to 
the initial singular state (H 2 — > oo). Meanwhile z = 
H = 0, which means in turn that Q = H5q = 0, i. e., 
the interaction vanishes. We have 



= 3iP 



1 



Accordingly, the autonomous system of ODE one obtains 
can be written in the following form: 



/V-^(l+2/ 2 -* 2 ) + ^ 



y 



-y\l + x l 



y 



2 I " V 3 

z' = 3z(z - + x 2 -y 2 ) 



H x 



(8) 



The 3D phase space where to look for critical points of 
(JSJ) can be defined as the following compact half-cylinder 
with unit height and unit radius: 

*3D = {{x,y,z)\0 < x 2 + y 2 < l,y > 0,0 < z < 1}. (9) 

We found five critical points P Ci : (x Ci , y Ci , z Ci ) £ ^3,0, 
which correspond to the five critical points found in 
21] for the non- interacting case. The stability proper- 
ties of these points, however, are modified by the non- 
gravitational interaction between the dark matter and 
the quintessence field. The main properties of the equi- 
librium points are listed below. 

1. Stiff matter-dominated solution(s), 

pf t : (±1,0,1), n v = 1. 



The quintessence field behaves like stiff matter, 
uj^p = 1, and the expansion is highly decelerated, 
q = 2. The eigenvalues of the Jacobian (also lin- 
earization) matrix, 



dx dx dx s 

dx dy dz 

dy' dy' dy' 

dx dy dz 

dz' dz' dz' 

dx dy dz 



o 



v0 



3t 







6, 



are: A x = 3, A 2 = 3 =p V^M, A 3 = 6. This 
means that, if /i < y/6, P s t is an unstable (source) 
equilibrium point, also called "past attractor" be- 
cause the orbits in the phase space approach to 
this point into the past. In case if [i > y/6, P^ t is a 
saddle critical point instead, which corresponds to 
absence of non-gravitational interactions (8q = 0). 
Meanwhile, P~ tl is always a past attractor. Since, 
z = 1, this solution can be associated either with 
the absence of additional non-gravitational inter- 
action (<5o = 0), or, with the initial singular state 
with diverging, H — > oo. 



2. CDM-dominated point, 

P m : (0,0,1), 0, 



1. 



This equilibrium point corresponds to decelerated 
expansion, q — 1/2, and the quintessence EOS pa- 
rameter, ujtp, is undefined. Since the eigenvalues 
of the corresponding linearization matrix, Ai = 
3/2, A2 = —3/2, A3 = 3, are of opposite sign, 
this means that this is a saddle critical point. A 
saddle in the phase space is " meta-stable" in the 
sense that it is unstable in some directions (x and 
z-directions in the present case), while being sta- 
ble in other directions in phase space (y-direction 
in this case). Orbits in which approach to the 
saddle point will evolve for some r-time in its neigh- 
borhood, until they are eventually repelled from it. 
This would mean that the saddle point corresponds 
to a transient stage of the cosmic evolution. In this 
case the matter-dominated phase is transient as it 
should be to account for the observed amount of 
cosmic structure. 

3. CDM/quintessence-scaling solution, 



p 

1 IT 




= 1 - 3/^ 2 

Q v 3/^ 2 



This equilibrium point exists whenever /1 2 > 3. It is 
associated with decelerating expansion (q = 1/2). 
The quintessence field mimics the behavior of pres- 
sureless matter (uj v = 0). The eigenvalues of the 
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corresponding Jacobian matrix are: 



A 1)2 = ^(-1 ± ^- 7 ),A 3 = 3. 

Hence, this is also a saddle point in ^30, so that 
the scaling solution describes a transient stage of 
the cosmic evolution. For /1 2 > 24/7, this is a spiral 
critical point. 

4. Quintessence-dominated solution, 

P q : Vl-M 2 /6, lj ^ % = 1, 

also known as potential/kinetic energy-scaling so- 
lution since the potential and the kinetic energies 
of the quintessence field scale as: 



V 



y 



1 - ^ 2 /6 



ip 2 /2 x 2 /i 2 /6 

This equilibrium point exists if /1 2 < 6. In this case, 

q = -1 + H 2 /2, u v = -l + M 2 A 

so that, the corresponding cosmological solution is 
accelerated whenever /1 2 < 2. Otherwise it is de- 
celerating. The eigenvalues of the Jacobian ma- 
trix corresponding to this point are: Ai = — 3 + 
/i 2 /2, A2 = — 3 + /j 2 , A 3 = [i 2 , so that, it is also a 
saddle critical point in ^30- 

Comparing our results with the results of reference 
[2l| . shows that, while there are not found new equi- 
librium points in the phase space in respect with the 
non-interacting case, the stability of the critical points 
P m - sc , and P q , has been indeed modified by the addi- 
tional (non-gravitational) interaction between CDM and 
quintessence components. Actually, while in the absence 
of interaction, depending on the value of the parameter 
/i 2 , the critical point P m _ sc can be, either a stable spiral 
{p 2 > 24/7), or a stable node (3 < \i 2 < 24/7), in the 
presence of constant interaction, Q = H5q, P m ~ sc is al- 
ways a saddle critical point. Similarly, in the absence of 
interaction the point P q is a stable node whenever /i 2 < 3 
and a saddle if 3 < [i 2 < 6. If we set-on the constant addi- 
tional interaction, Q = HSq, of non- gravitational origin, 
due to the fact that one of the eigenvalues, A3, is always 
positive, while the other, Ai, is always negative, P q is 
always a saddle point. 

Note that all of the critical points of the autonomous 
system of ODE ((5J, lay on top (z = 1) of the half-cylinder 
^3D ©: (x,y,l) € ^3D- !• e., there are not critical 
points associated with the additional non-gravitational 
interaction, Q = H5q. 



In the next subsections we shall look for modifications 
to the phase space dynamics of the present quintessence 
model which are induced by setting on an additional non- 
gravitational interaction between CDM and quintessence, 
in two generic ways found in the bibliography [12, HH, Ull , 
where S ^ const. 



B. 5 = 2,b 2 p m 

Here b 2 is a non- negative constant. In this case the x 
and y variables in Eq.© are enough and, besides: 

0- = ^ (l-x 2 -y 2 ) 

Hence, the autonomous system of ODE ([6]) can be written 
as: 



x' = + V 2 ~ x 2 ) + y 2 

3ft 2 

-2^(1-* -v) 



V> = \ V ( 1 + ^ ~ y2 ~ V 3^ 2 



(10) 



We shall seek for critical points of ([TO]) within the unit 
half-disc: 



*2D={(x,y)\0<x 2 +y 2 <l,y>l}. (11) 
Four critical points P Ci : (x Ci ,y Ci ) € Vt^D ar e found. 

1. Stiff matter solution(s): 

P±:(±1,0), u v = l,q = 2, 

which is dominated by the quintessence field since, 
f2„ = 1. The eigenvalues of the Jacobian matrix 
for Pf t are, Ai = 3(1 + b 2 ), and A^ = 3 T a/^A", 
respectively. Stability of these points is not af- 
fected in a significant way (with respect to the re- 
sults of |21|) by the non-gravitational interaction, 
Q = 3b 2 Hp m . 

2. CDM/quintessence scaling solution: 
b 2 ^ 




3(1 - 6 2 ) 2 



The point exists whenever, b 2 < 1, /1 2 > 3(1 — b 2 ) 
In this case CDM and quintessence scale as, 

n m _ i-(b 2 + 3(i-6 2 ) 2 / M 2 ) 

Q v b 2 +3(l-b 2 ) 2 /^ 2 ' 
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The quintessence EOS parameter, and the deceler- 
ation parameter are given by: 



P 2 b 2 



n 2 b 2 + 3{i-b 2 y 



q= 2 



2 : 



respectively. Note that, although both components 
of the cosmic fluid scale in a constant fraction, the 
quintessence does not exactly mimic the pressure- 
less CDM, since, lj v = 0, only in the absence of 
interaction (b 2 = 0). The eigenvalues of the lin- 
earization matrix are: 



Al,: 



"4(1 -b 2 )^ 



b 2 )(l 



b 2 )^ 



where we have defined, a := — 2/i 3 & 2 + 3(l- 
3b 2 )fi, and /3 := 4p 6 b 4 - 126 2 (1 - 6 2 )(5 - 
9(1 - 56 2 )(7 - 36 2 )^ 2 + 216(1 - b 2 f. Here we will 
not analyze in details the stability of this point but 
it is intuitive to see that the interaction does not 
modify in an appreciable way the results of (2lj | 
but for the fact that, if b 2 > 1/3, then q < 0, 
and the CDM/quintessence scaling solution may be 
correlated with accelerated expansion . 

3. Quintessence-dominated point (f2 v = 1), 



Pa 



(ji/VE, y/l - p 2 /&) ,U V = -1 + p 2 /3, 



which exists if p 2 < 6. The deceleration parameter, 
q = — 1 + /i 2 /2. The eigenvalues of the linearization 
matrix for this critical point are: Ai = — 3 + ^t 2 /2, 
A2 = —3(1 — b 2 ) + /j 2 . As seen, in this case, 
like in the former two ones, neither the stability, 
nor the existence conditions for the quintessence- 
dominated critical point are modified in any sensi- 
ble way by the additional (non-gravitational) inter- 
action, Q = 3b 2 Hp m . 

We can summarize this subsection by noting that, al- 
though the stability and existence conditions for the criti- 
cal points are not appreciably modified by the interaction 
term, Q = 3b 2 Hp m , the topology of the phase space it- 
self is indeed modified since one fixed point is missing 
with respect to the non-interacting case. Actually, it is 
easy to note that the CDM-dominated solution of the 
cosmological equations is not a fixed point of the equiv- 
alent autonomous system of ODE, as it is the case in 
the absence of interaction [2l|. The CDM/scaling point 
Pm-sc also differs from the corresponding one in the non- 
interacting case, since, i) the expansion may be acceler- 
ating as long as b 2 > 1/3, and ii) the quintessence EoS 
parameter lo v ^ 0, i. c., it does not track the CDM 
equation of state. This means that a mix of two different 



fluids: i) dust CDM, and ii) a scalar fluid with a different 
equation of state, co-exist during this stage of the cosmic 
evolution, in such a way that the ratio of their energy 
densities is a constant. 

A matter(CDM)-dominated point, or, in its absence, a 
matter/quintessence scaling point with io v = 0, is neces- 
sary to account for the observed amount of cosmic struc- 
ture. Hence, the addition of an interaction term of the 
form, Q — 3b 2 Hp m , to the quintessence field equations, 
transforms a suitable model of dark energy (as it is, in 
particular, the quintessence model driven by an exponen- 
tial potential) into a model which is ruled out by obser- 
vational data on structure formation. 



C. 5 = 3b 2 (p m +p v ) [2g|. 

Finally, to end up this 'illustrative' section, we shall 
briefly consider an interaction term which generalizes the 
former one, and has been studied in detail in several cos- 
mological contexts [26] . 

As before the phase space coincides with the unit half- 
disc, ^2D, of equation (JTTJ) . Since, Q/6H 3 = 36 2 /2, the 
autonomous system of ODE © simplifies to: 



1 t -t 2 2 

V = 2 V \ +x ~ y 



(12) 



Without going into details, one founds that, of the five 
equilibrium points found in the non-interacting case [5l[ , 
only three critical points remain: i,ii) the stiff matter 
solution(s), 




1 + VT+ 4b 2 







and, hi) the CDM/quintessence scaling point , 
p 



/ n 1 X He 

{Xc,yc), -p— = — o — , UJ U 



X c +Vc 



2 ■> 



1 1 

2 1 1 ' 



where, x c , and, y c , are the positive roots of the system 
of algebraic equations: 



1 + x 2 - y 2 - \l -fix = 0, 



x 2 (l + y 2 -x 2 )- x /-pxy 2 + b 2 = 0. 
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Hence, if consider an interaction term of the form 
Q = 36 H (p m + pu>), two important critical points are 
erased from the phase space by the non-gravitational in- 
teraction between CDM and DE: i) the (transient) CDM- 
dominated solution, which is necessary for the right for- 
mation of structure, and, ii) the quintessence-dominated 
solution which is important to get a stage of accelerated 
expansion which could be correlated with the present pe- 
riod of the cosmic history. 

In the next sections, on the grounds of a similar dy- 
namical systems analysis, we shall investigate the phase 
space dynamics of the so-called ghost dark energy mod- 
els, in order to seek for modifications to the correspond- 
ing dynamics originated from consideration of additional 
non-gravitational interactions. We shall be considering, 
in particular, the specific forms of the interaction dis- 
cussed in subsections III Al III Bt and III CI 



III. TIPS ON PHASE SPACE ANALYSIS 

There is an isomorphism between particular solu- 
tions of the cosmological field equations and points of 
the equivalent phase space. Take as an example, the 
quintessence cosmological model with exponential poten- 
tial studied in the former section, with additional non- 
gravitational interaction in the form of (subsection III Al) : 
Q = H5q. Consider the quintessence-dominated solution, 



P q : ( M /V6, y/l - ,i 2 /6, l) 



In this case, since, z = 1, then 
3H 2 



= 1 =► S = 0. 
Sti* + d 

Hence, this point will correspond to a particular so- 
lution of the quintessence field equations without non- 
gravitational interaction between the CDM and the DE. 
Besides, 



V6H V$ 



H 



a{<p) = aoe^, 



and, y 



VV/V3H = V 1 ~M 2 /6 which leads 



to 



y/2V /<p = \JQ — p? I p. The latter equation can be inte- 
grated in quadratures to give tp(t) = ipo+ln(t— £o) 2 i or, 
in terms of the scale factor: ait) = a (t — t ) 



a 



cto eg. The above solution for the original field vari- 
ables corresponds to the fixed point P q £ ~ the 
phase space defined in Eq.Q. This critical point exists 
if p 2 < 6 (see subsection III A[) . i. e., whenever the power 
2/p 2 > 1/3. The corresponding cosmological solution de- 
picts accelerated expansion whenever, p 2 < 2, i. e., if the 



power, 2//i 2 > 1. P q is a saddle equilibrium point. This 
means that there exists a non-empty set of initial con- 
ditions, which picks a congruence of neighboring orbits 
in ^3D, which approach to the quintessence-dominated 
point P q , then evolve for some r-time in the neighbor- 
hood of this point, until they are finally repelled from 
it. Going back to the standard cosmological "language" 
this would mean that the quintessence-dominated phase, 
depicted in spacetime by the above solution, can be, at 
most, a transient stage of the cosmic evolution. I. e., 
this cosmological solution can not decide the fate of the 
cosmic evolution. 

Studying the phase space associated with a given cos- 
mological model is, some times, more useful than looking 
for exact (and/or numeric) solutions of the cosmological 
field equations. Actually, while a given exact solution is 
necessary to evolve the model from given initial data to 
get concrete predictions to confront which with observa- 
tions, fixed points in the (equivalent) phase space tell us, 
which solutions in the (perhaps infinite) set of possible 
analytic solutions of the field equations, are generic, or 
" preferred" by the model. Generic or preferred solutions 
- corresponding to fixed points in the phase space - are 
structurally stable, while the remaining possible solutions 
are not stable in general. In this last sentence, under 
" structural stability" , we understand that the given crit- 
ical points are reached, either into the future, or into the 
past (in terms of the time-ordering phase space parame- 
ter t) quite independent of the initial conditions chosen, 
i. e., the qualitative behavior of the trajectories in the 
phase space is unaffected by small perturbations. In this 
regard, a unstable or source critical point can be alter- 
natively understood as a past attractor - a stable equi- 
librium point in reversed time direction - while a stable 
fixed point is a future attractor. 

We want to underline that, even if, according to the 
mentioned isomorphism, to every (physically meaningful) 
particular solution of the cosmological field equations of 
a given model, 3 it corresponds a given phase space point, 
there can be particular solutions of the field equations 
which can not be associated with fixed points in the phase 
space, i. e., which correspond to points in the phase space 
which are not critical. These solutions can be picked up 
only under very specific initial conditions on the phase 
space orbits, i. e., these will be structurally unstable so- 
lutions. In the subsection III Bl falso III Cl) . for instance, 
there are not found fixed points which could be asso- 
ciated with the matter(CDM)-dominated solution. This 



3 By " physically meaningful" we understand the fact that the given 
solution of the field equations corresponds to a point (a;, y, z) 



which belongs in the phase space if 



3D- 
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means that, even if this can be a particular exact solution 
of the cosmological field equations, it can not describe 
a matter-dominated stage of the cosmic expansion last- 
ing for enough time as to produce the observed amount 
of cosmic structure. Under the basis of this analysis, 
quintessential models with exponential potential, and, 
with additional non-gravitational interaction of the form 
Q = 'ib 2 Hp m (also, Q — 3b 2 H(p m + p v )) are ruled out. 

As we have just illustrated it in a simple quintessence 
model with an exponential self-interaction potential, a 
correct choice of the variables of the phase space is a very 
important point which is, most times, underestimated. 
One has to be able to choose variables which take val- 
ues within finite intervals, i. e., one has to be able to 
have the entire phase space enclosed in a finite region. 
This guarantees that all of the relevant fixed points will 
be "visible", i. e., none of them will go to infinity. For- 
tunately, most times, this choice is possible. Consider, 
for instance, an arbitrary phase space variable £ taking 
values in a semi- infinite interval: < £ < oo. After the 
following transformation of variable, 



i 



the above semi-infinite interval is transformed into a fi- 
nite interval < £ < 1. Using the above prescription 
one might transform an infinite interval ] — oo , oo [ into a 
finite one [— 1, 0[ 1J ]0, 1], by introducing two variables £ 
and £ to cover the entire phase space. 

In the remainder of this paper we shall replace the 
quintessence field by the so called QCD ghost dark en- 
ergy. We shall look for possible modifications to the GDE 
dynamics by an additional non-gravitational interaction, 
for the particular choices of the interacting term Q stud- 
ied in subsections III Al III Bl and III CI 



IV. PHASE SPACE DYNAMICS OF GHOST 
DARK ENERGY 

From now on we shall consider the dark energy com- 
ponent in the form of QCD (Veneziano) ghost [2h13||. We 
start with the simplest case when there is not additional 
non-gravitational interaction between the GDE and the 
CDM. For the flat FRW universe filled with GDE, CDM, 
and radiation, the corresponding Friedmann equation is 
written as 



3H = Pgde 



Pr, 



(13) 



where p m is the energy density of (pressureless) cold dark 
matter, p r is the energy density of radiation, and, Pgde is 
the GDE energy density assumed to be given by [4] , 



Pgde 



aH, 



(14) 



where, since we require non-negative p g d e > 0, and since 
we shall focus in expanding universes only (H > 0), then, 
the free parameter a > will be taken non-negative. 



A. Non-interacting Ghost Dark Energy 

The energy conservation equations for the different 
components are: 



p m + 3Hp m = 0, p r + 4Hp r ■■ 

Pgde + 3Hp gde {l + LOgde) = 0, 



0, 



(15) 



where uj gde is the GDE EOS parameter. The definition 
for the GDE energy density (fI5)l , together with Eq. (fT5j) , 
yield to the following relationship: 



Pgde H 

rj = ~FH = -3(1 + LOgde 

Hpgde H 2 



(16) 



It will be useful to have several quantities written 
in terms of the dimensionless parameter of GDE en- 
ergy density £l g de = Pgde I Hi 2 and of the dimension- 
less energy density parameter of the radiation compo- 
nent ft r = p r /3H 2 . In this regard, by taking the time 
derivative of the Friedmann equation (|13|) and, consider- 
ing equations ([TBI , and (IT5|) . one is left with 



3(1 + Ugde) 



3 3Qgde H~ 



H 
H 2 ' 



2 Qgde 

Hence, for the GDE state parameter one obtains 

^ 3(2 — Slgde)' 

while for the deceleration parameter: 



(17) 



(18) 



q = -1 + 3(1 + ugde) = 1 2 "y + ^ . (19) 

£ ^Lgde 

In order to put the cosmological equations in the form 
of an autonomous system of ODE, we choose appropriate 
phase space variables: 



X — Qgdej y — 



(20) 
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After this choice the following autonomous system of 
ODE can be obtained: 



3.t + y 



y' = -2y 



1 + x — y 



(21) 



where, as before, the tilde denotes derivative with respect 
to the variable r = In a(dr = Hdt). The phase space 
relevant to the present study is given by the following 
triangular (compact) region in (x, y)-plane: 



*2D = {{x ) y)\0<x<l, 0<y<l, 

0<x + y<l}. (22) 

In terms of the above variables, the Friedmann equation 
can be written in the form of the following constraint, 



Sl m = 1- x-y, (23) 
besides, the GDE EOS parameter can be written as 



3. GDE-dominated, de Sitter phase: 

p dS :(i,o), si m = o, si r = o, si gde = i. 

This late-time phase corresponds to an inflationary 
solution (q = — 1). The eigenvalues of the lineariza- 
tion matrix for Pds are: Ai = —4, A2 = —3, so that 
the solution is a future attractor. This equilibrium 
point mimics cosmological constant behavior since 

Ugde = — 1. 

According to the global structure of the phase space, 
orbits in i&2D converge towards the radiation-dominated 
stage into the r-past, depending on the initial condi- 
tions evolve for some time in the neighborhood of the 
saddle matter-dominated point, until they are repelled 
from this point to, finally, converge towards the GDE- 
dominated future (de Sitter attractor). This model cor- 
rectly describes the fundamental stages of the cosmic evo- 
lution comprised in the so-called 'present cosmological 
paradigm': i) a stage of radiation domination from which, 
ii) a period of matter-radiation equality and subsequent 
matter domination emerge, followed by iii) a late-time 
inflationary stage. 



Wgde — 



3(2 -x)' 



(24) 



Three equilibrium points P Ci : (x Ci , y Ci ) can be found in 
^2D, which correspond to different phases of the cosmic 
evolution: 



B. Interacting Ghost Dark Energy 

In the case when additional non-gravitational interac- 
tion between the GDE and the CDM is considered, the 
energy densities of GDE and CDM no longer satisfy in- 
dependent conservation laws. They obey, instead 



1. Radiation-dominated phase: 

p r :(o,i), a m = o, n r = i, si gde = o. 

This is a decelerating expansion solution (q = 1). 
The eigenvalues of the linearization matrix corre- 
sponding to this equilibrium point are: Ai = 2, 
A2 = 1, so that it is a unstable critical point (past 
attractor) in ^>2D ■ The GDE state parameter is 
Ugde = —1/3. 



2. CDM-dominated phase: 



P m : (0,0), Sl r , 



1, Sl r = 0, SI 



gde 



0. 



This phase of the cosmic evolution is characterized 
also by decelerated expansion (q = 1/2). The exis- 
tence of this solution is necessary for the formation 
of the observed amount of structure. The eigenval- 
ues of the corresponding linearization matrix are: 
Ai = —1, A2 = 3/2, so that it is a saddle critical 
point in ^2D ■ For the GDE state parameter we 
obtain uo g de = — 1/2. 



f) rn + 3Hp m = Q, p g de + 3Hp gd e(l + UJgde) = ~Q- (25) 

As before, for convenience, the interaction term can be 
written as Q = H8. In the present case equation (fT?| 
has to be replaced by the following: 



Hpgde 

In consequence one obtains 



Pgde 



3(1 + LUgde) = 

Hence, 

Ugde 



Qgde{3 — SSlgde + Sl r ) _ 2Sls 
Slgde{2 Slgde) 



Slgde{3 Sl r 



20, 



30g^ e (2 Slgde) 

1 — 2Sl g de + Sl r — Sis 



2-Sl 



gde 

H —3 4- 3Slgde — Sl r + Sl§ 

JP = 2 - Slgde ' 



(26) 



(27) 



(28) 
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FIG. 1: 3D "phase portrait" generated by given set of initial 
data for the autonomous system of ODE (|29[) . In the right- 
hand panel a view rotated 90° clock-wise about the z direction 
is shown. As shown in the left-hand panel, the phase space 
is bounded by the orthogonal planes x = 0, y = 0, z = 0, 
2 = 1, and by the oblique plane x + y = 1. It is seen the way 
the orbits of the ODE (|29[) depart from the global past at- 
tractor, point P r : (0, 1, 1)- radiation-dominated solution, are 
scattered from the saddle critical point P m : (0, 0, 1)- matter- 
dominated solution, and, depending on the initial data, these 
orbits approach to one of the points in the 1-dimensional man- 
ifold M s = {(x,0, 1/(4 - 3x)) : < x < 1} depicted by the 
dash-doted curve joining the critical points, Pas '■ (1, 0, 1) and 
P' m : (0,0, 1/4). The unphysical orbit drawn in gray color is 
originated by initial data leading to z-s smaller than 1/4 for 
a finite r-interval. 



Here, as in the former section, we will assume different 
forms of the interaction which are found in the bibliog- 
raphy [l2|, HH, H(| . We use same variables as in IIV Al 



S = So — const. 



Let us assume positive 5q > (energy transfer is from 
the GDE into the CDM). It is then appropriate to add a 
third variable (Eq.flJJ) 



3H 2 



3H 2 + S ' 



< z < 1, 



to the variables x, y in Ed.(|20[). The variable z takes 
decreasing values in the above interval. Since z = 1 cor- 
responds to the formal limit 5 — > 0, points of the phase 
space lying on the plane z = 1 are associated with non- 
interacting GDE solutions, unless the given fixed point 
can be associated with the initial big-bang singularity 
characterized by unbounded H 2 — > oo. The autonomous 
system of ODE obtained can be written as: 



z(4 - 3x + y) - 1 



y' = -2y 



z(2 - x) 
z(x — y) + 1 



z' = -2(1 - z) 



z(2 - x) 

z(4~3x + y) - 1 



2-x 



(29) 



The phase space in this case is the following 3D finite 
region (see FigfT]): 



*3D = {(x,y,z)\Q<x<l, 0<y<l, 

0<x + y<l, 0<z<l}, (30) 

where points lying on the plane z = 0, are not included 
since the ODE-s (|2"9"|) blow-up at z = 0. Here, as before, 
Q m — 1 — x — y, while 



Wgde 



x(3-y)z + 2(l -z) 



(y 



3x{2- 
2x)z + 2z - 



x)z 
1 



z{2 - x) 



(31) 



The equilibrium points, P Ci : [x Ci , y Ci , z Ci ), of the au- 
tonomous system of ODE ([29)) . together with their rele- 
vant properties, are listed below. 

1. (Non-interacting) Radiation-dominated solution, 
P r : (0,1,1). The cosmological parameters are: 
fl r = 1, q = 1 (ojgde is undefined). This is the past 
attractor in since the eigenvalues of the lin- 
earization matrix are all positive: Ai = 1, A2 = 2, 
and A3 = 4. 

2. (Non-interacting) CDM domination, P m : (0,0,1). 
In this case Q m = 1, q = 1/2 (uj g de is undefined as 
well) . As long as the eigenvalues of the linearization 
matrix for this critical point are: Ai = — 1, A2 = 
3/2, and A3 = 3, the solution is a saddle equilibrium 
point in ^30- 

3. 1-dimensional manifold, 



{(z, 0,1/(4 -3a;)) : < x < 1} . 



whose points can be associated with interacting 
CDM/GDE-scaling solutions whenever, < x < 1, 
since the CDM and the GDE scale as, n m /£l g d e — 
(1 - x)/x. The point (1,0, 1) VL gde = 1, which 
also belongs in the manifold, can be associated, 
instead, with a non-interacting GDE-dominated 
phase, while the equilibrium point, (0,0,1/4) => 
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Q m — 1, corresponds to an another phase of mat- 
ter domination (see below). 

For points in the manifold M s , the deceleration pa- 
rameter q — — 1 (accelerated expansion), while the 
GDE EOS parameter is given by, u> g de = —1/%, so 
that, as x decreases from 1 to 0, ui g de unbound- 
edly decreases from —1 to -co. Hence, this mani- 
fold can be associated with phantom-like solutions. 
Equilibrium points in M s are non-hyperbolic crit- 
ical points since one of the eigenvalues of the lin- 
earization (also Jacobian) matrix 



_ 3x 
2-x 
x 
2-x 
x(4-3xf 
2-x 





-4 





18(1-2:) 
(4-3x) 2 (2-x) 

6(l-x) 
(4-3x) 2 (2-x) 
6(l-g) 
2-x 



is vanishing: Ai = 0, A2 = —4, and A3 = —3. Nev- 
ertheless, the manifold itself is stable in the sense 
that every orbit in the phase space end ups in M s 
no matter which initial conditions were given (see 
the figure [1]). This manifold joints the the non- 
interacting GDE-dominated solution (ft g d e = 1) 
Pds '■ (1,0,1), with the interacting CDM domina- 
tion critical point (fi TO = 1) P' m : (0, 0, 1/4). As P' m 
is approached, the GDE state parameter blows-up, 



lim u> gde 

(x,j/,z)-K0,0,l/4) 



—00, 



so that this phase of the cosmic evolution can be 
associated with extreme phantom behavior. 

The phase portrait for this case (S — 5q) is depicted 
by the figure [TJ As seen, phantom-like behavior is the 
generic fate of the cosmic evolution, while a de Sitter 
universe corresponds to only one (boundary) point in the 
manifold M s : the non-interacting GDE-dominated solu- 
tion Pds- Although the introduction of an interaction 
term of the form Q = H5o does not modify the stabil- 
ity properties of the fixed points as compared with the 
non-interacting case QIV A[) . nevertheless, the phase space 
topology is indeed modified since the isolated fixed point 
corresponding to GDE domination in the non-interacting 
case is replaced by the stable manifold M s . Points in 
this manifold can be associated with phantom-like solu- 
tions of the cosmological equations without the need for 
phantom-like matter. As a consequence, for this particu- 
lar form of the interaction term Q, since the phase space 



4 Recall that, with the course of expansion, since H 2 decreases 
from infinitely large values in the past, into vanishingly small 
values in the future, in the manifold M s , z = 3H 2 / (3H 2 + 8q) 
decreases from 1 to 1/4, so that x = (4z — l)/3z takes decreasing 
values in the interval [1,0]. 



dynamics is richer, the additional non-gravitational in- 
teraction makes the model more attractive. 

In the following cases (subsections IIV B2] and IIV B~3~|) . 
the phase space coincides with the 2D phase plane cor- 
responding to the triangular region defined in equation 

(El. 



2. S = 3b 2 p m . 

As before b 2 is a positive constant. Here SI5 = 36 2 (1 — 
x — y). The autonomous system of ODE obtained for this 
particular case is given as follows: 



3(l-6 2 )(l-x) + (l + 3fe 2 )y 



y' = -2y 



2-x 

(l + 3b 2 )(l-y) + (l-3b 2 )x 
2-x 



(32) 



For the GDE EOS parameter ui g d e and the deceleration 
parameter q — —(1 + H/H 2 ), one obtains 



J gde 



x(3-6b 2 -y) + 6b 2 {l-y) 
~ 3x(2-x) ' 

1 - 36 2 - (2 - ib 2 )x + (1 + 3b 2 )y 



(33) 



2-x 

The critical points for this case are: 

1. Radiation-dominated solution 

P r :(0,l), Q r = l, q=l. 

It is the past attractor in (Ai = 2, A 2 = 

1 + 3b 2 ). The expansion is always decelerated in 
this regime. The GDE-EOS parameter uj g d e is un- 
defined. 

2. CDM-dominated solution 

P m :(0,0), fi m = l, q = (1 - 36 2 )/2. 

Here u) g d e is also undefined. This phase of the cos- 
mic expansion is accelerated whenever, b 2 > 1/3. 
For this case the eigenvalues of the Jacobian ma- 
trix are, Ai = 3(1 - 6 2 )/2, and, A 2 = -1 - 36 2 , 
so that, it is a saddle equilibrium point if, b 2 < 1. 
Otherwise, for b 2 > 1, it is the future attractor in 
the phase space, ^20- 

3. GDE-dominated (de Sitter) solution 

PdS ■ (1,0), Vl g de = 1, Ugde = -1, Q= -1- 
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Since the eigenvalues of the Jacobian matrix for this 
point are Ai = —3(1 — b 2 ) and X 2 = —4, it is the 
future attractor whenever b 2 < 1, while for b 2 > 1 
it is a saddle critical point instead. 

Addition of the non-gravitational interaction, Q = 
3b 2 p m , between CDM and GDE, does not modify the 
topology of the phase space in any significant way. The 
only effect of the additional interaction is to slightly mod- 
ify the stability properties of the already existing fixed 
points in the non-interacting case. However, a very inter- 
esting thing happened: the CDM-dominated phase can 
be associated with accelerated expansion if b 2 > 1/3. 
Besides, if b 2 > 1, the stage of matter dominance not 
only supports accelerated expansion, but it is the future 
attractor in the phase space. In this case the GDE- 
dominated solution represents a transient stage of the 
cosmic evolution. 5 



3. 5 = 3b 2 (p m + p gde )- 

In this case Cl$ = 3b 2 (I — y). We obtain the following 
autonomous system of ODE: 



3(l-6 2 ) + (36 2 + l)y-3x 



y 



-21/ 



2-x 
(3b 2 + 1)(1 - y) + x 



2-x 



(34) 



As in the former case, here the phase space is the same 
as in Ea.(|22p. The magnitudes of relevance can be writ- 
ten in terms of the phase space variables (fi m = 1 — x— y): 



Ugdt 



q = 



x(3-y) + 6b 2 (l-y) 
3x(2-x) 
l-3b 2 - 2x+(l + 3b 2 )y 



(35) 



The following equilibrium points are found: 

1. Non-interacting, radiation-dominated phase 

P r : (0,1), n r = 1, a, = 0, 

while ujgde is undefined. Since the eigenvalues of 
the linearization matrix, Ai = 2, A2 = 1 + 36 2 , are 
both positive, this is a unstable critical point (a 
past attractor), which corresponds to a decelerating 
expansion solution (q = 1). 



2. CDM-dominated phase 

P m :(0,0), O m = l, q = (l-3b 2 )/2. 

As before the GDE-EOS parameter oj g de is unde- 
fined in this case. The eigenvalues of the lineariza- 
tion matrix corresponding to this critical point are: 



Ai 



-1 - 3b 2 , A 2 = -3/2(6 2 - 1), so that this 



solution is a saddle critical point if b < 1, while, 



for 6 > 1, it is the future attractor in, ^2D- 
b 2 > 1/3 the solution is inflationary (q < 0). 



For 



3. CDM/GDE-scaling solution P m - sc : (1 - 6 2 ,0). 
This equilibrium point exists whenever b 2 < 1 (for 
b 2 = 1 this point coincides with P m above) . In this 
case the energy densities of DM and of GDE scale 
as H, m /ftgde — b 2 /(1 — b 2 ). It is a late-time infla- 
tionary solution (q = —1). The eigenvalues of the 
linearization matrix for P m _ sc are: Ai = —4, A2 = 
3 — 6/(1 + fe 2 ), so that, since b 2 < 1, this solution 
is the late-time (future) attractor. The GDE-EOS 
parameter uj gde = l/(6 2 — 1) —00 < uj gde < — 1. 
I. e., in this regime the GDE always behaves like 
phantom matter. Since, according to the observa- 
tional evidence, at present 



< 



< 1 => 



gde 



l-b 2 



< 1, 



then b 2 < 1/2. 



Consideration of an interaction term of the form Q = 
3b 2 H(p m + Pgde) does actually modify the topology of 
the phase space ([22]) by replacing the GDE-dominated 
(de Sitter) solution in the non-interacting case by a 
CDM/GDE scaling solution. Hence, the additional non- 
gravitational interaction between the cold dark matter 
and the ghost dark energy, in this case, smooths out the 
coincidence problem. Actually, since the scaling solu- 
tion, whenever it exists, is the (global) future attractor, 
then, independent on the initial conditions chosen, scal- 
ing between dark matter and dark energy always dic- 
tates the late-time cosmological behavior. Besides, at 
the CDM/GDE scaling point the expansion is acceler- 
ated since q = — 1. 



V. GDE PLUS A DECAYING A 

Here we consider a further generalization of the GDE 
model explored so far (Eq.([H|)). which is given by [6] (see 
also 0) 



5 It would be interesting to show that the GDE-dominated critical 
point can be associated with early-time inflation for b 2 > 1. 



Pgde = aH + /3H 2 , 



(36) 
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where, in order to warrant positivity of energy, the con- 
stant parameter /3 > should be a non-negative quan- 
tity. The second term in the above expression is inspired 
in models of decaying cosmological constant 19]. One 
may expect that the subleading term H 2 in (|3"S|) might 
play a crucial role in the early evolution of the universe 
(20l | . It was shown that considering this term can give 
better agreement with observational data compared to 
the usual GDE [61. 



A. Non-interacting model 

The cosmological field equations are the same as in 
IIV Al (equations fp~3|) . (|T5|) ). Here we use same variables 
x and y defined in (I2U1) . Notice that the relationship (|3"6")l 
can always be written in the following alternative form: 



H 



3x- p. 



(37) 



It is seen from (|3T|) that, since, 

1 / a\ 
*=g(/J + ^),0<*<l, 

then H is bounded from below: 

0<P + %<3^H>H = ^. 

This entails that the de Sitter fixed point can describe the 
late-time cosmological dynamics in this model, while the 
decaying cosmological constant can drive the early-times 
dynamics since, near of the initial singularity 

P 

H ->• oo x —,=$■ a^O => pgde^pH*. 
o 

Given that the free parameter a > is non-negative 
and, since we are interested only in expanding solutions 
H>0, then, from ((37J it follows that 

— = 3x- 8>0 x > -. 
H H ~ ~ 3 

In consequence, the phase space of this model will corre- 
spond to the compact triangular region <~ $2D- 



%Z D = {(x,y)\0/3<x<l, 0<y<l, 

0<x+y< 1}. (38) 

Due to (|36l equation ([16) is replaced by the following 
equation: 



Pgde (3x + 0\ H 

- ' JH = -3(l + w ffde ). (39) 



Hp g d 



3x I H 2 



The corresponding autonomous system of ODE is now 
given by, 



03 - 3z) 
"3(1 



3{x-l)-y 



6-P-3x 
h x - y) - 2/3 
6-P-3x 



(40) 



As before Q m — 1 — x — y, while the expressions for the 
GDE-EOS and the deceleration parameters, are modified 



Ugde 

q = 



1 



(3x + p)[3(x-l)-y] 



3x(6 - f3-3x) 
3 + P - 6x + 3y 
6-P-3x 



The following fixed points are found in ^^d - 

1. Radiation/GDE-scaling solution 

P r _ sc : (/9/3,l-0/3), q = l. 

The GDE EOS parameter uj g d e = 1/3 whenever 
P 7^ 0. For P = it is undefined instead. This point 
exists whenever p < 3. Given that the eigenvalues 
of the corresponding Jacobian matrix are positive: 
Ai = 2, A2 = 1, it is a global past attractor (source 
point) in ^^d- The ghost dark energy tracks ra- 
diation with their dimensionless energy density pa- 
rameters scaling as 

n r 3-0 



Qgde P 

2. CDM/GDE-scaling solution 

P m - sc :{P/3,0), uj g de = 0, q=l/2, 

which exists only for < 3. This case is similar to 
the former one: ghost dark energy tracks cold dark 
matter and 



gde 



3-/3 



but, this time, it is a saddle critical point since the 
eigenvalues of the linearization matrix for this point 
are of opposite sign: \i = 3/2, A 2 = — 1. 

3. de Sitter (GDE-dominated) phase 

Pds ■ (1,0), uj g de = -1, q= -1- 

This is an inflationary solution. The ghost dark 
energy behaves like a cosmological constant. This 
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solution is a global (future) attractor in the phase 
space $2Di since both eigenvalues of the corre- 
sponding Jacobian matrix: Ai = —3 and A2 = —4, 
are negative. This means that Pds is the end-point 
of any phase space orbit and, hence, it depicts al- 
ways the late-time cosmological dynamics. This is 
consistent with the analysis of the expression (|37|) 
at the beginning of this subsection, according to 
which the Hubble parameter H, is always bounded 
from below (H > Hq). 

Take a look at ([37]) . it is then clear that, for points 
where x = (3/3 (fixed points P r - sc and P m - sc ), the ghost 
dark energy behaves like a decaying cosmological con- 
stant pgde = 0H 2 , while for the GDE-dominated de Sit- 
ter solution Pds, since x = 1, the GDE behaves like a 
standard (truly constant) cosmological constant: 



H = H 



3-/3 

Pgde = aH + /3Hq = 



»- ^ =A. 



(3-/9)5 



B. Interacting Model 

We keep using the same phase space variables, x = 
Qgde, y — fir- In this case the conservation equations 
(fT5|) are replaced by (P2"5j) (the conservation equation for 
radiation is unchanged). The relation (f39|) . distinctive of 
ghost dark energy, transforms into the following relation- 
ship, 



Pgde 
H Pgde 

where Q$ 



3x + /3 



3x 

S/3H 2 . Besides, 



m = - 3(1 



- Ugde ) — j 



(41) 



1 i? _ 3(x - 1) - y + il s 
3H 2 ~ 6 - (3 - 3x 



Also, 



Wgde 



■1 - 



(3z + 0)[3(a:-l)-tf]+6fi 4 



(42) 



(43) 



3x(6 - j8 - 3x) 
A very general form of the interaction term can be 

S = Sq+ 3b 2 (p m + Pgde), 

which comprises in a single expression the interaction 
terms of III A\ III Bt and III Cl However, as in the former 
sections, it is convenient to study separately the different 
particular cases. 



1. 5 = Sn 



In the present case the already defined variables 



X — Qgdei V — f^r; Z 



3H 2 



3H 2 + 6 ' 



are a good set of phase space variables in order to write 
the cosmological equations (p~3|) , (|25|) , plus the conserva- 
tion equation for the radiation component, in the form 
of an autonomous system of ODE. In terms of these vari- 
ables the dimensionless energy density associated with 
the interaction £1$ = fl$ = Sq/3H 2 , can be written as 



1-2 



(44) 



while the GDE-EOS and the deceleration parameters can 
be written as 



Ugde 



! _ (3x + /3)[3(a - 1) - y]z + 6(1 - z) 
3x(6 — (3 — 3x)z 



and 



-1-3 



[3(a? -l)-y]z+l-z 
(6-/3- 3x)z 



respectively. The autonomous system of ODE for this 
case is the following: 



x' = (/3- 3x) j 



-2y 



[3(x-l)-y]z+l- z 
(6 - /3 - 3x)z 
3(1 + x - y)z + 1 - z - 2/3z 



(6-/3- 3x)z 



while the phase coincides with the 3D compact region 



$t D ={(x,y,z)\0/3<x<l, < y < 1, 

<x + y< 1, < z < 1}. (46) 

The following critical points/manifolds can be found 
in the phase space ^l D : 

1. (Non-interacting) radiation/GDE-scaling point, 

P r _ sc : 08/3, 1 - 0/3, 1) , q = 1, 

where the dimensionless energy densities scale as, 

fl r _ 3 - /3 

figde P 
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Whenever (3^0, the GDE EOS parameter uj gde = 
1/3 (if P — it is undefined). This fixed point 
exists if (3 < 3. This is either a non- interacting 
point, or it is associated with the initial big-bang 
singularity instead (H 2 — > oo is unbounded) since 
z = 1 =>■ (5 /3H 2 ) -> 0. The eigenvalues of 
the linearization matrix for this point are: Ai = 1, 
A2 = 2, and, A3 = 1/3, so that P r - sc is a global 
past attractor in ^!^ D . 

2. (Non- interacting) CDM/GDE-scaling fixed point, 

P m -sc : 08/3, 0, 1) , uj gde = 0, q = 1/2, 

which exists whenever, /3 < 3. The CDM and GDE 
dimensionless energy density parameters scale in a 
constant fraction, 

3 p 



(a) (Interacting) CDM/GDE-scaling fixed point, 



and GDE tracks CDM since, u) g de = w m = 0. Given 
that, the eigenvalues of the Jacobian matrix for, 
P-m-sci are °f opposite sign: Ai = 3, A 2 = 3/2, 
A3 = —1, this is a saddle critical point in ^l D . 

3. 1-dimensional compact manifold, 



where 



x,0, 



1 



4- 3x 







< x < 1 > , 



Jgde 



-1/x, q = — 1. Besides, 

Qrn 1 X 



Since (the real part of) one of the eigenvalues of the 
Jacobian matrix for points in this manifold, P Si € 
M s , vanishes: 



Ai — 0, A2 — —3, A3 — — 



4(3 - p) 



3 - /3 + 3(1 - x) ' 



then, the P Si , are non-hyperbolic critical points. 
This means that no final judgment can be made 
about their stability as isolated points. We want 
to underline, however, that the manifold M s itself 
is stable in the sense that, no matter which initial 
conditions to choose, every orbit in the phase space 
will, necessarily, end-up in a given point belonging 
in M s . 

For, /3/3 < x < 1, fixed points, P Si £ M s , corre- 
spond to (inflationary) phantom-like (uj g d e < — 1), 
interacting CDM/GDE-scaling solutions. 

The end-points of this manifold are: 



/ ,: 1-0 1 

m-sc ■ [ 3,0, A _ f3 



Ugde — -"g, q = — I, 



for which, fl m /il g d e = (3 — /3)//3. Unlike 
the (non-interacting) matter-scaling solution, 
P m - Sc , since, in this case, q = — 1, this fixed 
point is associated with inflationary expan- 
sion. Besides, even if the cold dark matter and 
the ghost dark energy scale in a constant frac- 
tion, due to the additional non-gravitational 
interaction, the GDE does not track the CDM: 

3 

U g de = 7^ w m = 0- 

(b) (Non-interacting) de Sitter solution, 
Pds = (1, 0, 1) , uj g de = -1, q= -1. 

Since, £l g d e — 1, this is a GDE-dominated so- 
lution. In this case, according to (f3"T)) . 



H = H, 







3-/3 



There exists another 1-dimensional critical manifold, 

6 \V H, 7-f3-3x) 3" " 3 
for which, 

Q m = [3 - 2 - X, UJgde 



6-/3 
3x 



, q=-l. 



However, for this manifold to be physically meaningful, 
it is required that, i) f3 < 3, and, ii) (6-/3)/3 < l. 6 Both 
constraints, together taken lead to, j3 — 3. Therefore, 
the only physically meaningful fixed point in M* is the 
de Sitter solution, PdS: already included in the manifold 



2. S = 3b' 2 p m 

The phase space for this and the following cases of 
interacting CDM and GDE is the two-dimensional tri- 
angular region, defined in Eq. (|38l) . Since, in the 
present case, 5 = Zb 2 p m , we have, 



The latter constraint comes from requiring, x < 1. 
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so that, 



Ugde 



-1-3 



fi 4 = 3b 2 {l-x-y), 



(3x + P)[3(x - 1) - y] 
3x(6 - P - 3x) 

66 2 (1 — x — y) 
x(6 — /3 — 3x) ' 

3(x-l) -y + 36 2 (l - as -y) 



(47) 



6 - P - 3x 



We are led to the following autonomous system of 
ODE: 



x' = (P - 3x) 
-6y 



3(x~l)-y + 3b 2 {l 



y) 



y 



6-/3 

l + x-y + 3b 2 {\- 



3x 



y) 



3 



6 - p - 3x 



(48) 



This time, basically the same three fixed points of sub- 
section [VA] can be found in the phase space how- 
ever, their stability properties, among others, are slightly 
modified by the interaction: 

1. Radiation/GDE-scaling solution, 



03/3,1-^9/3), q=l, 



where 



Qgde 



3-/3 
P ' 



For p 7^ 0, ojgde = 1/3- This fixed point exists if 
P < 3. The eigenvalues of the Jacobian matrix are 
both positive: Ai = 2, and, A 2 = 1 + 3b 2 , so that, 
as in the non-interacting case, the point, P r _ sc , is 
a past attractor in the phase space. 

2. CDM/GDE-scaling phase, 

P m _ sc : (P/3, 0) , iu gde = -3b 2 /p, q=(l- 36 2 )/2, 

which, as with the former critical point, exists 
whenever p < 3. As in the non-interacting case, 
the energy densities involved scale in a constant 
fraction: 



n 



gde 



3-/3 
~P~' 



However, this time the ghost dark energy does not 
track the cold dark matter, since, in general (b 2 ^ 



0, P ^ 0), to gde ^ uj m = 0. Besides, for b 2 > 1/3, 
this solution depicts accelerated expansion, while, 
for b 2 > P/3, the ghost dark energy behaves like a 
phantom field (without the unwanted features of a 
phantom field, of course). 

The stability is also modified by the interaction in 
respect to the non-interacting case, since the eigen- 
values of the corresponding linearization matrix: 
Ai = 3(1 - b 2 )/2, A 2 = -3(1 + b 2 ). Hence, when- 
ever, b 2 > 1, the CDM/GDE-scaling (inflationary) 
solution is the future attractor in the phase space. 

3. de Sitter solution, 



PdS ■ (1,0), ttgde = 1, Ugde = ~1, <7 



-1. 



In this case the additional non-gravitational inter- 
action, S = 3b 2 p m , affects only the stability prop- 
erties, since, Ai = —4, and A 2 = —3(1 — b 2 ). It 
is seen that, for b 2 > 1, the de Sitter solution is a 
saddle critical point in while, for b 2 < 1, it is 

the future attractor instead. 



3. 8 = 3b 2 (p m + pgde) 

As before we use variables x = £l g d e , and y — Q r , tak- 
ing values in the compact triangular region, Q38[) . In 
this case the dimensionless 'energy density' of the interac- 
tion, fl s = 3& 2 (l-y). For the GDE-EOS and deceleration 
parameters we get: 



Ugde 



(3x + P)[3(x - 1) - y] 
3x(6 - P - 3x) 

6b 2 (l-y) 



-1 



x(6 — P — 3x) 
3{x-l)-y + 3b 2 {l-y) 



Q-P-3x 

We obtain the following autonomous system of ODE: 



x = (P - 3x) 



3(x-l)-y + 3b 2 (l-y) 



Q-P-3x 



-6y 



l + x-y + 3b 2 {l-y) 
6 - P - 3x 



2n 



(49) 



The non-gravitational interaction, in this case, appre- 
ciably modifies the topology of the phase space with re- 
spect to the non-interacting case. The equilibrium points 
found in, ^ 2 _d> are the following: 
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1. CDM/GDE-scaling solution, P m _ sc : 09/3,0), 
% 2 3-/3- % 2 



0(3 - 0) 



2(3-/3) 



This fixed point exists whenever ft < 3. The eigen- 
values of the linearization matrix are: 

A - 1 %2 A - 3 (l 3b2 \ 

so that P m _ sc is stable if 3(1 - b 2 ) < ft < 3. Oth- 
erwise it is a saddle. 

2. Interacting CDM/GDE-scaling solution, 



p L~ sc : (1 - fo 2 ,0), w grfe 



1 -fe 2 



which exists if b 2 < 1. The eigenvalues of the lin- 
earization matrix are: 



Ai = -4, A 2 



3(3-0 -3b 2 ) 
3-0 + 36 2 ' 



This critical point is stable whenever, either ft < 
3(1 - b 2 ), or ft > 3(1 + 6 2 ). If 3(1 - b 2 ) < ft < 
3(1 + b 2 ) it is a saddle point. 



There is another critical point, 

'ft 3~ft + 9b 2 



P 



r — rn — sc 



3' 3(1 + 3& 2 ) J ' 



3b 2 

but, in this case, il m = — prggs < 0, which is not physi- 
cal, so this point does not belong in ^^d- 



interaction term Q may lead to CDM/quintessence scal- 
ing critical point which is correlated with accelerated ex- 
pansion if, either Q = 3b 2 Hp m , or Q = 3b 2 H(p m + p v ). 
For the case when Q = H5q the matter-scaling solution 
is associated with decelerated expansion as it happens in 
the non-interacting case. 

At the same time the non-minimal coupling of CDM 
to the quintessence field leads to unwanted effects. De- 
pending on the kind of interaction considered the mod- 
ifications to the dynamics might be either: i) minimal 
if Q = HSq since only the stability is modified but not 
the asymptotic structure of the phase space, or ii) sub- 
stantial if Q = 3b 2 Hp m or Q = 3b 2 H(p m + p v ): the 
asymptotic structure is modified by the interaction since 
one or several critical points are erased. In the two lat- 
ter cases the CDM-dominated point is erased, while the 
CDM/quintessence scaling solution is transformed from 
stable into a saddle critical point in the phase space, 
i. e., matter-scaling results in a transient stage of the 
cosmic evolution. Besides, at the matter-scaling point, 
even if the energy densities of the dark matter and of 
the ghost dark energy are related in a constant fraction 
Qcdm/Qgde = const, the quintessence field does not track 
the CDM equation of state parameter: lo v ^ u m — 0. 
This means that the resulting model may not take ac- 
count for the observed amount of cosmic structure. The 
following conclusion may be read off from the former re- 
sult: since the additional non-gravitational interaction 
Q modifies the phase space topology in the previously 
discussed way, a suitable model of dark energy is trans- 
formed by the interaction into a model that may be ruled 
out by the observations on structure formation. 



A. Ghost Dark Energy 



VI. DISCUSSION 

Consideration of non-minimal coupling between the 
matter degrees of freedom and the scalar field, or, al- 
ternatively, of additional non-gravitational interactions 
between matter and the dark energy, has been repeat- 
edly invoked to get late-time matter-scaling attractors 
which allow accelerated expansion [22-24]. A scalar field 
is expected to couple explicitly to ordinary matter unless 
some special symmetry prevents or suppresses the cou- 
pling [22, 28]. For the case studied in section|TT]of this pa- 
per, i. e., scalar (quintessence) field with an exponential 
potential, even in the minimal coupling situation (stud- 
ied in Ref . [2l1] ) , there exists a matter-scaling late-time 
attractor, however it is a decelerated solution in contra- 
diction with the present cosmological paradigm. As we 
have shown here, consideration of the non-gravitational 



The situation is a bit different if the quintessence field 
is replaced by a QCD ghost. In this case, in general, 
the interaction makes the GDE model more attractive 
and better suited to accommodate the present cosmolog- 
ical paradigm without the coincidence problem. In the 
simplest situation when the additional non-gravitational 
interaction is not considered, there are 3 equilibrium 
points in the equivalent phase space of the GDE model. 
For the model p g d e = olH these are: i) the radiation- 
dominated point P r : (0, 1) - the past attractor, ii) CDM- 
dominated critical point P m : (0, 0) - saddle equilib- 
rium point, and iii) de Sitter (GDE-dominated) phase 
Pds ■ (1,0) - the future attractor. Meanwhile, for the 
model pgde — olH + ftH 2 the critical points are: i) 
the radiation-scaling point P r - S c ■ (0/3,1 — 0/3), ex- 
ists if < 3, the GDE mimics radiation Lu g d e = 1/3, 
Vl r /Slgde = (3 — ft)/ ft - the past attractor, ii) CDM- 
scaling point P m - sc ■ (0/3,0), also exists if ft < 3, the 
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ghost field tracks cold dark matter io g d e = <^cdm = 0, 
^m/Qgde — (3 — (3) / (3 - a saddle point, and iii) GDE- 
dominated (de Sitter) point Pds '■ (1,0) - the future at- 
tractor. As seen, the role of the term oc H 2 is to re- 
place radiation domination and matter domination by 
radiation scaling and matter scaling respectively. This is 
expected since the QCD ghost, in this case, is not vanish- 
ing at any stage of the cosmic evolution: at early stages 
the term oc (3H 2 dominates, while at late times it is aH 
which dominates the cosmic dynamics. The topology of 
the phase space is also modified by the quadratic term 
since there is bifurcation at the critical value (3 — 3: for 
f3 < 3 the 3 critical points can be found, but for (3 > 3 
only the de Sitter point survives. 

Let us discuss now what the role of the additional 
non-gravitational interaction is. The constant interac- 
tion S = 5q modifies mainly the GDE-dominated de Sitter 
point in both cases, 6 — 3b 2 p m and S = 3b 2 (p m + p g de), by 
replacing it by a one-dimensional manifold M s , which is 
characterized by scaling of CDM and GDE: Sl m /f2 sc je = 
(1 — x)/x. This manifold is correlated with accelerated 
expansion (q = —1) and joints a matter-dominated (ac- 
celerated) phase (x = 0), with a GDE dominated one 
(x = 1). Since uj g d e = —1/x, points in the critical man- 
ifold are associated with phantom QCD GDE. Although 
critical points in M s represent scaling of the dark matter 
and of the ghost dark energy, which represent acceler- 
ated expansion, the coincidence problem is not allevi- 
ated by the interaction Q = HSq. Actually, even if every 
phase space orbit, no matter what the initial conditions 
are, ends up at a point in M s , the given constant ratio 
Qm/Qgde is highly dependent on the initial conditions. 

If consider the GDE model with p g d e = aH, an in- 
teraction term of the form Q = 3b 2 Hp m modifies the 
properties of the matter dominance point P m : (0,1), 
Q m — 1, which is associated with decelerated expansion 
(q = 1/2) in the uncoupled case, not only by allowing 
a range in parameter space b 2 > 1 where matter domi- 
nance is the late-time attractor, but also by allowing a 
range b 2 > 1/3 where accelerated expansion may be asso- 
ciated to dark matter domination. For b 2 > 1 the GDE- 
dominated de Sitter solution Pds ■ (1,0) (u g d e = — 1, 
q = —1) is a saddle critical point in the phase space 
(it is the future attractor otherwise). In the case when 
the interaction term Q = 3b 2 H(p m + p g de) is assumed, 
the modified matter-dominated critical point shares the 
same properties with the case of interaction 5 — 3b 2 p rn 
(future attractor for b 2 > 1, associated with acceler- 
ated expansion if b 2 > 1/3). In addition the de Sitter 
phase in the uncoupled case is replaced by matter-scaling 
Pds Pm-sc : (1 - b 2 , 0), ft n /n gde = b 2 /(l - b 2 ), which 
is associated with accelerated expansion (q = —1), and 
exists if b 2 < 1. It is the late-time attractor whenever 
it exists. Hence, in this long as b 2 < 1, the co- 



incidence problem is actually resolved by the additional 
non-gravitational interaction. 

For the QCD GDE model with p gde = aH + (3H 2 , 
in the case of interaction Q — 3b 2 p m , the stability of 
the matter-scaling point P m _ sc , Cl m /Q g d e = (3 — (3)/ (3, 
is modified for b 2 > 1, in which case P m _ sc can be the 
late-time attractor. Besides, for b 2 > 1/3 this point is 
correlated with accelerated expansion and a resolution 
of the coincidence problem is also achieved. Notice that, 
since uj g d e — —3b 2 / (3 ^ 0, the ghost dark energy does not 
mimic dark matter (y> c dm = 0). It is also noticeable also 
that the matter-scaling point exists whenever j3 < 3. The 
situation is not as attractive if consider the interaction 
term Q — 3b 2 {p m + p g d e )- In this case only CDM-scaling 
solutions survive: the radiation-scaling point is erased, 
while the de Sitter point is replaced by a new phase of 
matter scaling P d s -t P* n - sc ■ (1 - b 2 ,0). We are left 
with two different phases of dark matter/ghost dark en- 
ergy scaling behavior: i) P m _ sc : ((3/3,0) which exists if 
(3 < 3 and is the future attractor if 3(1 — fo 2 ) < (3 < 3, and 
ii) P^n_ sc : (1 — b 2 , 0) which exists if b 2 < 1 and is the late- 
time attractor if either (3 < 3(1 - b 2 ), or (3 > 3(1 + b 2 ). 
P m - sc is correlated with accelerated expansion whenever 
3(1 - 3b 2 ) < (3 < 3, while P* n ^ sc depicts always accel- 
erated expansion since q = — 1. Even if the coincidence 
is resolved in this case by the non-minimal coupling, the 
resulting model might be ruled out by the observational 
evidence on structure formation since there are not crit- 
ical points that can be correlated neither with radiation 
nor with dark matter dominance. 



VII. CONCLUSIONS 

We have performed a detailed study of the phase space 
of the QCD ghost dark energy models where p g d e = aH 
and p g de = aH + (3H 2 respectively. Both cases with and 
without additional non-gravitational interaction between 
cold dark matter and the Veneziano ghost field have been 
scrutinized. Special attention has been paid to the cor- 
rect choice of the phase space variables which allow for 
bounded and compact spaces in each case. As a result all 
of the possible critical points of physical relevance were 
found. It has been shown that the non-minimal coupling 
introduces more feasible changes in the asymptotic struc- 
ture of the equivalent phase space of QCD GDE models 
than for quintessence models with an exponential poten- 
tial. 

In this paper we have neglected the contribution from 
the dynamics of the trapping horizon as it is done in most 
works on GDE 

[Ml- This particular case, strictly 
speaking, corresponds to dark energy dominance so that, 
in general, the latter contribution to the dynamics has to 
be taken under consideration. We defer investigation of 
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this more realistic case for further work. 
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